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Measurements of the two dimensional angular correlation of the electron-positron annihilation 
radiation have been done in the past to detect the momentum spin density and the Fermi surface. 
We point out that the momentum spin density and the Fermi Surface of ferromagnetic metals can be 
revealed within great detail owing to the large cancellation of the electron-positron matrix elements 
which in paramagnetic multiatomic systems plague the interpretation of the experiments. We prove 
our conjecture by calculating the momentum spin density and the Fermi surface of the half metal 
Cr02, who has received large attention due to its possible applications as spintronics material. 

PACS numbers: 71.18.+y, 71.27.+a, 71.60.+Z, 78.70.Bj 

Keywords: density functional theory, positron annihilation, 2D- AGAR, Cr02, Fermi surface 



To a great extent, the Fermi surface (FS) can regarded 
as the defining property of a metal. It is an ubiquitous 
concept which appears in numberless works devoted to 
the study of the electronic structure of systems which, un- 
der some particular circumstance, seem to show metallic 
behaviour. In several cases, the initial task of establishing 
metallic behavior is not easy since different probes can 
yield contrasting answers. Clearly, an experimental in- 
vestigation of the FS implies the observability of the elec- 
trons sitting in the partially filled energy bands, which, 
in turn, requires interaction of the conduction electron 
with the experimental probe. 

In the case of a measurement of the two dimensional 
angular correlation of electron positron annihilation ra- 
diation (2D- AGAR), the FS is revealed through disconti- 
nuities (breaks) in the electron positron momentum den- 
sity, p^^(p) at points = {kpi + G^), where G is a 
reciprocal lattice vector and are the reduced Fermi 
wave vectors in the first BZ. While the locations of the 
FS breaks are faithfully preserved 2], the resulting sin- 
gle particle electron momentum density, (p) , is severely 
modulated by the non uniform positron (spatial) density 
and the electron-positron Coulomb interaction. In this 
letter we show that in the special case of a ferromagnetic 
sample, the probe modulation is strongly suppressed, al- 
lowing a faithful representation of the spin polarized mo- 
mentum density and unambiguous interpretation of the 
data. The FS breaks are reinforced by the Lock-Crisp- 
West (LCW) transformation 0, consisting of folding the 
momentum distribution p^^(p) back onto the first BZ by 
translation over the appropriate vectors G. The result of 
the summation (denoted as LCW density) is 0, 

n ^ 

(1) 

Here and (j) denote the electron and positron wave 
function, respectively, E^? is the Fermi level and ek,n is 



the energy eigenvalue of the electron from band n with 
Bloch wavevector k. The enhancement factor |^ 5'(r), 
describes the enhancement of the electronic density at 
the positron location due to the Coulomb force. If the 
k, n dependence of the overlap integral in eq.Qis negligi- 
ble, p^^^ (k) is reduced to the occupancy^ i.e. the num- 
ber of occupied bands per k point. The FS manifolds 
are the loci of the breaks of the occupancy. A signifi- 
cant drawback of 2D- AC AR is that all the outer electrons 
overlap (to some extent) with the positron probe. There- 
fore, when the orbital character of one full valence band 
changes noticeably in the BZ, the related LCW density 
acquires a k-dependence which is superimposed to the 
changes in the occupancy due to the FS breaks of the con- 
duction bands. For example, our recent work on UGas 
has detected a noticeable change in the p character 
of three valence bands located 1-2 eV below which 
did obscure greatly the visibility of the FS. Further no- 
table examples are high-Tc superconductors where strong 
positron wave function effects prevented the observation 
of the critical FS sheets linked to the Cu-Oxide planes, 
which are responsible of superconductivity 0, • 

Obviously, the elimination of the contribution of all 
the filled bands from the LCW summation would increase 
greatly the visibility of the FS whenever the overlap inte- 
gral of eq. Qis /c(valence)-dependent. This favorable sit- 
uation is indeed realized in the measurement of the spin 
polarized bands of ferromagnetic metals. In this case, 
the employment of 2D-ACAR experiments 0, Q 
hinges on two facts: i) the intrinsic polarization Pg+ for 
positrons produced during (5 decay (for the ^^Na source 
Pe+, averaged over angle and velocity, is about 36%). 
ii) the annihilation selection rule which requires that the 
positron may undergo 27 annihilation only if the spins of 
the annihilating pair form a singlet state. 

Therefore, by performing 2D-ACAR measurements 
with the magnetic substance polarized (by an external 
magnetic field) in directions respectively parallel and an- 
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tiparallel to the average positron polarization (which is 
unchanged upon reversal of the magnetic field) will then 
yield an observable imbalance of 27 annihilations with 
respect to the majority or minority spins. The exper- 
imental spectra taken when the positron polarization is 
parallel (antiparallel) to the polarizing magnetic field will 
be 



P7ar.anU,aM = 1/2(1 ± Pe^W^ {v) + iV^^(p)/2]- 

+1/2(1 T Pe^)[N^{v) + N^^(V)I% 



(2) 



where refers to the contribution from the spin 

degenerate (i.e. non magnetic) bands. Subtraction of 
p^^^ from p^antipar yields the net momentum spin density 



Pllinip)' 



P^pm(P) =Pe+[p7arip) 



P antipar (P) 



(3) 



Equations [21 and 01 are equally applicable to the electron 
positron momentum density, p^^(p), or to the LCW den- 
sity, yielding for the latter case, the LCW net spin density 

I LCW spin 

(k). 

The main point of this work is the disappearance of 
and all the positron wave function effects related 
to the bands in question. As obvious as it may appear, 
this interesting result has, to our knowledge, never been 
pointed out and should prove to be very useful in future 
studies of ferromagnetic metals. 

To investigate our conjecture we have chosen to exam- 
ine Cr02 which is predicted to be a half metal by a va riety 
of ab initio band structure calculations 0, 0, Q, Il5j . 
Due to this result Cr02 has attracted much attention in 
the field of spintronics for its potential use as injector of 
a highly (nominally 100%) spin polarized current in the 
future field effect spin transistor. A puzzling result of the 
calculations is that different recipes to approximate the 
exchange correlation potential, U^cc (local spin density 
approximation fLSDA) Il6l and Generalized Gradient 
Approximation (GGA) [I3) lead to unusually noticeable 
differences in the FS topology A further interest- 

ing result is that the standard LSDA and GGA seem to 
explain ultraviolet photoemission spectroscopy (UPES) 
experiments [isj better than the LDA+[/ method, 0, 
which usually is more suited to treat strong electron cor- 
relations. 

In Refs. 0, we have presented a method to calculate 
directly the LCW density via eq. Q As a base for our 
calculation we used the scalar relativistic full-potential 
linearized augmented plane waves method (FP-LAPW) 
implemented in the WIEN2k pack age j2oj . Compared 
to similar works of other authors |j,l2l|, our scheme, in- 
cluding electron-positron correlation effects IBj , being full 
potential and compatible with any option of WIEN2k, 
which include spin-orbit (s.o.) interaction, LDA+/7 and 
orbital polarization, is more suited to the study of nar- 



row bands and electron correlations. 

The extension of the scheme to spin polarized calcula- 
tions is relatively straightforward. The same steps used 
in the paramagnetic calculations are done for spin up 
and spin down electrons separately. In the absence of 
s.o. interaction the result consists of pure up spin and 
down spin bands. After the inclusion of s.o. each Bloch 
state becomes a combination of spin up and down compo- 
nents. For each band crossing E^? there is now an up and 
down spin submanifold, both having the same occupancy 
(and therefore the same Fermi surface). The amplitudes 
of these submanifolds are, however, different. Owing to 
parity and angular momentum conservation, the positron 
essentially projects the electronic states over the spin up 
or down manifolds. 

The crystal structure of Cr02 is simple tetragonal (Ru- 
tile type, space group P42 /mnm, ^^) with two formula 
units per unit cell. Details of the calculations, performed 
adopting \]xc according to GGA and LDA in the param- 
agnetic and ferromagnetic phases, are reported elsewhere 
[23. The calculated majority spin bands are in excellent 
agreement with those reported by Mazin et al The 
main finding of the calculation is the half metallic be- 
haviour, with a large energy gap (:^2eV) in the minority 
spin bands. The conduction bands have a very strong Cr 
d character, which is ascribed to states of t2g symmetry. 
It is worth noting that the Cr-d orbital character of all 
the conduction bands is rather constant along each of the 
high symmetry directions. 

As mentioned above, the topology of the holelike FS 
sheet, shown in Fig. mi,II)b depends noticeably on the 
kind of Ua,c (LSDA,GGA). Whereas in the GGA case this 




FIG. 1: The two FS sheets (a,b) of FM Cr02, shown in 
half BZ (Online color: a:orange, electronlike; b:blue, holelike) 
produced with U^c from LSDA (I) and GGA (II) E 113. 
Here and in the next figures capital letters label the high 
symmetry points of the BZ. 

FS has a simple, pillow-like shape (Fig.p/II-b), the corre- 
sponding FS manifold resulting from LSDA (Fig. ^/I-b) 
is a pseudosphere connected to a toroidal structure with 
rhombic shape. On the other hand, the electronlike, F 
centred structure, yielded by the two calculations is very 
similar. Since Cr02 is a compensated metal (a necessary 
condition of any half metal) , the electronlike and holelike 
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sheets have equal volume. The Fermi volumes resulting 
from LSDA and GGA differ slightly, corresponding to 
12.8% of the BZ and 10.7% of the BZ for LSDA and 
GGA, respectively. 
The predicted response of a 2D- AGAR experiment ex- 




FIG. 2: (Online color.) FLAPW calculations for Cr02 in the 
(010) plane of the tetragonal BZ: (a) up spin LOW density; 

(b) down spin LOW density; (c) (a) + (b); (d) (a)-(b). The 
black contour denote the intersections of the two main FS 
sheets (see fig^ with the (010) plane. Note that the orizontal 
direction is along [001] 

periment can be elucidated by a slice of the calculated 
LGW density (adopting LSDA) in the (010) plane (GGA 
yields little difference in this plane; recall Fig In 
Fig. [2 panels (a) and (b) refer to the majority and mi- 
nority spins LGW densities, respectively, whereas panels 

(c) and (d) denote sum and difference of panels (b) and 
(b). These panels would be obtained by an experiment 
with 100% spin polarized positrons impinging a sample 
where the polarizing magnetic field is parallel or antipar- 
allel to the positron polarization. In panel (a) the breaks 
pertaining to the electronlike and holelike FS sheets ap- 
pear rather clearly (compare the changes of gray scale 
(or color) with the contour line marking the FS breaks), 
in spite of a strong modulation caused by a non uniform 
positron density. This modulation is the only feature 
appearing in panel (b), consistent with a negligible spin 
up-down mixing due to s.o. (WIEN2k yields 99.97% elec- 
tron polarization at Ei?). In panel (c) we have summed 
panels (a) and (b) to notice that positron wavefunction 
effects are reinforced at the expense of the FS signatures. 

The very appealing message, and the main point of 
this paper, is however provided by panel (d), where most 
of the positron wave function effects cancel out and the 
almost intact topology of the FS is restored. 

As anticipated, the difference eliminates the contribu- 
tion from the non spin polarized bands and any positron 
wave function effect to be ascribed to them. The can- 
cellation of positron effects caused by all the non mag- 
netic bands, which is obviously effective also for ferro- 
magnetic insulators , is clearly of great importance in 




FIG. 3: Spin polarized FSs seen by 2D-ACAR experiments 
in Cr02, shown in half BZ (theoretical prediction). The FSs 
are identified by the two isodensity surfaces (a,b) extracted 
from p^^ cr,,>.(k) as described in the text. (Online color: 
a:orange, electronlike; b:blue, holelike). The pertaining LCW 
densities were produced with U^c from LSDA (left, I) and 
GGA (right, II) 



systems with several atoms in the unit cell and, conse- 
quently, several bands near to Ep. A similar effect can 
be inferred from Fig. 9 of Ref.[ll|a), referring to the half 
metal NiMnSb. In Gr02 the visibility of the FS is further 
enhanced by the weakness of the positron modulation of 
the conduction bands. This result is consistent with the 
constancy of the d orbital character in the BZ for all the 
conduction bands, noted above. 

The similarity of panel [2td) with the theoretical oc- 
cupancy (whose breaks denote the FS) makes feasible 
algorithms aimed at extracting the FS from isodensity 
surfaces of the LGW density. This task would be impos- 
sible for the data set shown in panell^^a), where, owing to 
positron modulation, the amplitude of the LGW density 
nearby the Fermi break of a single sheet changes notice- 
ably in the BZ. Conversely, p^^^ spinO^) Presented in 
panel Eljd) allows to separate very clearly the two con- 
duction bands and apply the method presented in Ref. 
ji^ which identifies a multi-sheet FS in terms of iso- 
density surfaces selected at the loci of maximal ampli- 
tude variation of the LGW density. It is worth noticing 
that in the half metals the analysis of p^^^ spinO^) 
the search of the FS manifolds is particularly suited. In 
fact, since the subtraction of the insulating LGW (down) 
spin density from the conducting LGW (up) spin den- 
sity preserves the proper sign in the jumps of the oc- 
cupancy due to the FS breaks, maxima of p^^^ spinO^) 
will correspond to maxima of the occupancy and vicev- 
ersa. Fig |31 shows the resulting isodensity surfaces of 
l-^Tcw spinO^)- "^^^ similarity with the true FS shown 
in Fig. Qis striking. The cancellation of positron wave 
function effects shows that 2D- AGAR has the power to 
investigate the FS predictions of LSDA and GGA, shown 
in Fig. n The difference in the FS is well revealed in 
the (110) plane. To establish to what extent feasible ex- 
periments can accomplish this task, in Fig. ^ we have 
simulated in full the output of the experiment, accord- 
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FIG. 4: (Online color.) FLAPW calculations for Cr02 in the 
(110) plane of the tetragonal BZ. (a) LSDA up spin occupancy 
convoluted with the experimental resolution; (b) LSDA sim- 
ulated experimental difference spectra (inclusive of statistical 
noise; see text), ^^^^^^^^(k); (c) same as (a) for GGA; (d) 
same as (b) for GGA. 



ing to eq. [3 The steps adopted to produce Fig. 21 are, 
in succession: i) Calculation (adopting X^xc from GGA 
and LSDA) of p^^paA^) and p"^^^ antipaA^) accord- 
ing to the mixing shown in eq. [21 employing the realis- 
tic positron polarization Pg+ = 0.35; ii) Convolution of 
P?cwpar{^) and p'^^^antipar(k) with the experimental 
resolution R [25|; iii) Perturbation of the two simulated 
densities with statistical noise. It is assumed that the 
LCW density was reconstructed from 5 projections, each 
collecting 5x10^ coincidence counts, and that the spec- 
tra were symmetrized in the 3D reconstruction procedure 
0; iv) Construction of p^/^^ ,p,„ (k) . 

Figure ^ shows the electronic occupancy (reflecting 
the FS topology) convoluted with the experimental res- 
olution [panels (a),(c)] compared to the difference spec- 
tra resulting from steps i)-iv) [panels (b),(d)] for LSDA 
and GGA, respectively. Interestingly, the bulging of the 
LSDA holelike FS at 2± (0.27,0.20) a.u. [panel (a)], absent 
in the GGA holelike FS [panel (c)] is well reproduced by 
the difference /)^^^g^^^(k) [panel (b)] and not present in 
panel (d). Obviously, to reveal the FS subtleties at stake, 
arrays with large number of counts are required, particu- 
larly when differences between experimental spectra are 
analyzed. With this caveat, we predict that a 2D- AGAR 
experiment should decide over this issue. 

In conclusion, we have pointed out the power of the 
magnetic LGW procedure to eliminate large part of the 
modulation of the LGW density due to the non uniform 
positron density and reveal the FS in great detail. We 
have proved that this indeed happens by applying our 
scheme to calculate the LGW density of the half metal 
Gr02. We have simulated in full the output of an experi- 



ment performed with partial polarization of the positron, 
dictating the amount of data collection required to have 
appropriate signal to noise ratios. Note that the ability 
to include s.o. effect can be employed to corroborate 2D- 
AGAR measurements aimed at determining the polar- 
ization of the conduction electrons at the Fermi surface. 
This type of information could be of critical importance 
for the design of novel materials to be employed in the 
implementation of the future spin-based FET. 
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